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Abstract 

We investigate the group field theory formulation of the EPRL/FK spin foam models. 
These models aim at a dynamical, i.e. non-topological formulation of 4D quantum 
gravity. We introduce a saddle point method for general group field theory amplitudes 
and compare it with existing results, in particular for a second order correction to 
the EPRL/FK propagator. 
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1 Introduction 

Group field theories (GFTs) |T] are quantum field theories over group manifolds and can 
be also viewed as higher rank tensor field theories [21 E] which generalize matrix models. 
They provide one of the most promising frameworks for a background invariant theory of 
quantum gravity in which one sums both over topologies and geometries. Indeed, each 
Feynman graph of a D dimensional GFT can be dually associated with a discrete space- 
time via a specific triangulation and gluing rules given by the covariance and vertices of the 
theory. The functional integral formalism defines a weighted sum over triangulations with 
each weight (amplitude) related to a sum over geometries via a spin foam formalism (3j (see 
[3 E] for results on power counting and non-perturbative resummation of such models). 

Spin foams are the Feynman amplitudes of GFT. But GFT in addition specifies the 
class of graphs that should be summed, together with their combinatoric factors. This 
stems from Wick contractions rules, hence (perturbative) GFT requires to distinguish the 
non-quadratic part (interaction) from the quadratic part (propagator) in the field action. 

The simplest group field theories correspond to quantization of the BF models, hence 
to topological versions of gravity. Recently new spin foam rules have been proposed for the 
quantization of full fledged 4D gravity [TJ, El El [Hi] ■ These models stem from an improved 
analysis of the Plebanski simplicity constraints. The corresponding so-called EPRL or 
FK models are neither of the BF nor of the Barret-Crane type. They mix the left and 
right part of 5*0(4) ~ SU(2) x SU(2) in a new way which gives a central role to the 
Immirzi parameter. These new theories could be called dynamical since their propagators, 
combining two non-commuting projectors, have non-trivial spectrum. 

Preliminary studies of the asymptotic large spin (also called "ultraspin" ) regime have 
been performed for the EPRL/FK amplitudes of the "self-energy" and the "starfish" graphs 
(see Figures [5] and [3]) [TTJ. These results are a first step towards a study of renormalizability 
of such theories. 

In [12] a linearized approximation has been devised to investigate the ultraspin limit 
of BF spin foam amplitudes (see also [T3] ). This approximation captures the correct power 
counting of some graphs, such as type 1 graphs in the Boulatov model [H] , but it typically 
overestimates more general graphs. 

In this paper we push further the group field theory approach to the EPRL/FK mod- 
els, first introduced in [9], and perform another step towards the general investigation of 
their renormalizability. We use a coherent state representation of the EPRL/FK propaga- 
tor, as in [5], while other representations are exhibited mainly for comparison with other 
approaches. We introduce a general saddle point approximation, as in [13], which repro- 
duces correctly the approximation [12J to the power counting of BF amplitudes for sim- 
ply connected graphs and, for non- degenerate configurations, the EPRL/FK "self-energy 
graph" power counting of [11]. We discuss also the case of degenerate configurations, not 
studied in [TT] . 

The plan is as follows: the next section is devoted to a review of the BF and EPRL/FK 
group field models in a field theoretical spirit. The following section presents the stationary 
phase method. Finally we remark that the sign of the self-energy graph points towards a 
singularity in the effective propagator of the EPRL/FK model, which could signal a phase 
transition. For completeness we included useful formulas and normalization conventions 



in the appendices. 

2 Implementation of GFT 

In GFT, the field arguments live on products of Lie groups. Feynman amplitudes are 
spin-foams, i.e. two-complexes with vertices, stranded lines (also called propagators) and 
faces (that is, closed circuits of strands). 

2.1 Fields 

Since GFT represents a quantum theory of space-time itself, the usual spin-statistic the- 
orem may not apply. In this paper we consider only Bosonic statistics. However other 
choices have been considered [TB]. We also work with an Euclidean signature. 

The number of strands in the GFT lines encodes the space-time dimension D. The 
natural group associated to such a D dimensional GFT is [SO(D)] D , hence a field is 
a function on [SO(D)] D . We don't assume any symmetry under permutations of the 
arguments. 

2.2 Vertices 

In the spin-foam literature the term "vertex" usually refers to the vertex together with 
the square roots of its dressing propagators. This terminology is not the standard one in 
quantum field theory. Further confusion often stems from the fact that in BF theory the 
propagator is a projector hence is equal to its square, and also to its square root! 

To clarify the situation, let us return to ordinary field theory. In that case also 
the definition of the vertex could be considered ambiguous since one can dress it with a 
more or less arbitrary fraction of the propagator. What fixes this ambiguity is the usual 
requirement that vertices in field theory should obey a certain locality property in direct 
space. This allows to distinguish them from their dressing (half)-propagators, which are 
non-local operators. 

Since GFT is non-local on the group we cannot transpose directly this rule. To 
properly distinguish the vertex from the propagator we propose to use an extended notion 
of locality adapted to the GFT case, which we call simpliciality [T7] . 

For consistency reasons every vertex in GFT is required to have a total degree in the 
fields ensuring parity of the number of strands. In odd dimensions this restricts the degree 
of the vertex to be even. Hence we propose the following definition: 

Definition 2.1 A vertex joining 2p strands is called simplicial if it has for kernel in direct 
group space a product of p delta functions matching strand arguments, so that each delta 
function joins two strands in two different half-lines. 

The usual vertex for D— dimensional GFT is a <p D+l simplicial vertex in which the 
faces are glued in the pattern of a D-dimensional simplex. For instance the ordinary 




Figure 1: A simplicial vertex of a 4— dimensional GFT. We have chosen here a particular 
matching and orientation for each of the strands. 

Boulatov vertex in 3 dimensions is simplicial (with p = 6) as it writes 
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satisfying to our definition. But remark that the "pillow term" [5] 
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£3)0(£3, 34, 35)0(35, 34, 3e)0(3e, 32, 3i)- (2.3) 



is also simplicial in D = 3. Also in any dimension D there are infinitely many higher than 
order D + 1 simplicial vertices according to our definition. 



2.3 Propagators 

We consider only field theories in which the propagator C is Hermitian. It can be consid- 
ered either as an Hermitian operator —¥ C0 acting on fields or as its Hermitian kernel 

C(gi,...,g D ;g'v-,g' D )' 



[C0](3i,---,3d) = / dg' 1 ...dg' D C(g 1 ,...,g D ;g' 1 , 



!,-••, sf D ). (2-4) 



The corresponding normalized Gaussian measure of covariance C is noted d/ic- Hence 



C(3i,...,3d;3i,...,3z)) = / 0(3i,...,3d)0(3i,...,3d)^c. 



(2.5) 



2.4 Graphs 



Graphs are generated by gluing together propagators and vertices, according to Wick 
contractions, hence to Feynman rules. 
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Figure 2: A stranded propagator with particular orientation; two strands have rjtf = 1 and 
the other two have r\nf — — 1. 

Definition 2.2 A stranded graph is called regular if it has no tadpoles (hence any line £ 
joins two distinct vertices) and no tadfaces (hence each face f goes at most once through 
any line of the graph). 

It is convenient to introduce orientations on both lines and faces of stranded graphs. 
Regular oriented graphs are natural since they are conveniently described by two matrices 

• the ordinary incidence matrix e v j. which has value +1 if the edge £ enters the vertex 
v, —1 if the edge £ exits vertex v and otherwise. Hence Y2 V \ e v,t\ — 2 for each t. 

• the incidence matrix rjij between faces and edges, which has value +1 if the face / 
goes through edge £ in the same direction, —1 if the face / goes through edge £ in 
the opposite direction and otherwise. Hence ^\ \r]£j\ = D for each £ (see Fig. [2|. 

These orientations are useful to write down the integrand of the Feynman amplitudes. 
However the integrals, that is the spin-foam amplitudes themselves, do not depend on these 
orientations, at least for the class of theories considered in this paper. 

From now on we consider only amplitudes for regular graphs. This is for convenience, 
as generalization to any graph of our formulas is possible. It has been argued that GFT 
should in fact be restricted to colored models, which generate only regular graphs [TS1 HEJ 
[19]. Remark that every colorable stranded graph is regular, but the converse is not true; 
colorable graphs in particular have all their faces of even length, hence the "starfish" graph 
of Figure [3] with ten faces of length 3, although regular, is not colorable. 

2.5 The BF theory 

2.5.1 Propagator in direct space 

By direct space we mean the representation which uses the group elements. 

In the case of the BF theory the propagatoiQ here noted P, is just the projection on 
gauge invariant fields: 



. ..,d D ) = dh (f)(gi h,...,g D h), (2.6) 

where the integral is performed over the group SO(D) with respect to its Haar measure. 
Let us remark that P 2 = P so that the only eigenvalues are and 1 (which means that the 

1 Beware that this propagator is called C in [6| [20] . 




Figure 3: The "starfish" graph, quantum correction to the vertex. The dashed lines rep- 
resent the edges (they do not correspond to strands). Each edge contains 4 strands, there 
are 40 such strands, forming 10 closed faces and 10 open faces. We have shown 4 faces: 3 
closed and one open which take into account all 4 strands of one particular edge (upper 
left). 



BF theory has no dynamics). The operator P is Hermitian with kernel 



W{gi,-,9D;g[,-g' D ) = / dh JJ<5(^/i(^) _1 ). 



(2.7) 



2.5.2 Amplitudes in direct space 

Suppose that we choose an arbitrary orientation of the lines and faces of a graph Q (which 
for simplicity has no external legs). 



Combining together the vertices (2.1) (or the generalization to dimension D) and 



the propagators (2.7) of the graph, the integration over all g variables can be explicitly 



performed, leading to the direct space representation of the BF Feynman amplitude as an 
integral over line variables h: 



J e&L g feF g v 7 

where Lg, Fg are the set of lines and faces of Q, respectively. The oriented product 
Yliefh 7111 means that the product of the variables he has to be taken in the cyclic ordering 
corresponding to the face orientation (starting anywhere on the cycle). 

As announced, the amplitude ( 2.8[ ) neither depends on the arbitrary orientation of 
the lines, nor on those of the faces. A pedestrian way to see this is to exploit carefully the 
parity of the 5 function and of the Haar measure under h — > h^ 1 . Beware that formula 



(2.8) may be formal as this amplitude can be infinite for many graphs. 



2.5.3 Amplitudes in the angular momentum basis 



The angular momentum basis uses the irreducible unitary representation spaces V J ' of 
dimension dj = 2j + 1. In this space there is a standard decomposition of unity 



h = s ^\^ m ){h 



m\ 



(2.9) 



where \j,m),m G [— j, j] is the usual orthonormal basis in V 3 . 

In dimension 3, using the Peter- Weyl theorem we can transform (2.8) to get the 
representation: 



where the Ponzano-Regge vertex is the 6j symbol (cfr appendix A) corresponding to the 
6 face indices meeting at the vertex. 

In dimension 4 for simplicity we work with the covering group SU(2) x SU(2) of SO(A) 
and we decompose the group elements as g = (g + , g_), with g± e SU{2). Moreover we write 
j = (j + ,j_) for the eigenvalues of the angular momentum J in each SU(2) component. 
The Peter- Weyl decomposition leads to the similar angular momentum representation: 

4 = I]E d 3f+ d 3f „l[{l5j + }{l5j-}. (2.11) 

/ if+Jf- v 

where {15j} are the 15j symbols (see for example pj] for the definition and normalization 
conventions). 



2.5.4 Coherent states 

Consider ^'"Kg), the matrix element of the group element g in the representation j, 
computed between the states (j,m\ and \j, k). We have 



l i = d i y)|i, m)(j, ro'| / dgR^ m J (g)R U) J m '(g) = d 3 [ 

JSU(2) JSU(2) 

where we have introduced the notation: 



\j,g)(j,g\, (2.12) 



(2.13) 



The states \j, g) are the coherent states |21j . and the last expression in (2.12) is a decom- 
position of the identity in terms of these coherent states. 



Let us recall that the decomposition of the identity (2.12) can be further simplified 
and taken over the coset G/H, G = SU(2), H = U(l): 



h = dj / dn\j,n){j,n\ 

Jg/h=S9. 



(2.14) 



with \j, n) = g n \j,j) and g n denned in (A. 15). We suppress the domain of integration G/H 
in what follows. 

The states \j,n) form a generating set in V^ sometimes called "overcomplete basis". 

Let us now turn to the coherent states of the group SU(2) x SU(2). In fact one has 
four possible such coherent states which are given by acting with the same group element 
(#+) 9-) £ SU(2) x SU(2) on either of the following states: 

\j,j)®\j,j), \j,j)®\j,-j), \j,-j)®\j,j), \j,-j)®\j,-j)- (2-15) 

Note that these four states can be obtained from one another by the action of an SO (A) 
group element. However, if one considers only the action of the diagonal SU(2) subgroup 
of elements of the form (g, g) then there are two inequivalent states that cannot be related 
by such a transformation. 

2.5.5 The BF propagator and amplitudes using coherent states 

To prepare for the EPRL/FK propagator we rewrite the BF propagator inserting the 
coherent state decomposition of identity on each strand. Let us consider SU(2) BF first. 
Since P 2 = P we can introduce two distinct SU(2) gauge-averaging variables, u and v at 
both ends of the propagator, instead of the single variable h (e. g. u on the side where 
e Vj £ = —1 and v on the side where e Vy g = +1). Between these two variables we insert the 



partition of unity (2.14). This does not modify the propagator. Working out the algebra, 
we find 

r 4 

/) = / d^nE d ^ T %' («0/W)-v' ii i/) > ( 2 - 16 ) 

J /=i if 

with gf,g'f,u and v elements of SU{2). The index / labels the four strands of the propa- 
gator, which belong to four different faces (since we consider only regular graphs). 

To write down the amplitudes we need to introduce some notations. There are now 
group variables 2|Lg|, U£ and vg, and D\Lg\ variables n. The amplitude is again factorized 
over faces: 

Ag= f \{ du^vz \{ A f . (2.17) 

To write down Af, let us number the vertices and lines in the (anti)-cyclic order along a 
face / of length p as £i, V\ • ■ ■ £ p , v p , with by definition £ p+1 = £ 1 . We have then 

A f = E < +1 / n dn ^f < j' n ^K:L h Z a + tl b'. n ^j >' ( 2 - 18 ) 



where 



he a ,v a = v ia if e Va/a = +1 

h a ,v a = u ia if e Va/a = -1. (2.19) 



2.5.6 The D = 4BF case 

In D = 4 we work with SU(2) x SU(2), the covering group of SO(A); we have the similar 
system \j + ,n + ) ® \j_,nJ) of coherent states and the partition of unity on the space Vj, 
with j = (j+,jJ) 



lj_ = lj = d j+ dj_ / dn + dn^\j + ,n + ) ® |j\_,n_)0+,n+IO'-> n -l- (2.20) 

The gauge- averaging variables, u = («+,«_) and v = (v + ,v_) at both ends of the prop- 
agator are now elements of SU(2) x SU(2). Between these two variables we insert the 



partition of unity (2.20) and we find 

i 



g')=fdudvf[ J2 ^f^f-^v^v^Amia'fr'v- 1 !^) (2.21) 

/=lif_i_,.7>_ 



with gf,g'f,u and v elements of SU(2) x SU(2), and we have formulas similar to (2.17)- 



(2.18) for the amplitudes. 



2.6 The EPRL/FK GFT 

The EPRL/FK model introduces a modification of the propagator of the BF model, while 



the vertex remains the same. The EPFL/FK propagator has a structure similar to (2.21) 
but with replacement of lj by a non-trivial projector. We notice at this point that since 
this projector does not commute with P, it is not possible to recombine u and v in a single 
gauge averaging variable h. 

It implements in two steps the Plebanski constraints with a non-trivial value of the 



Immirzi parameter 7. Starting from the (2.21) expression of the BF propagator in the 



coherent states representation, the first step adds the constraint j+/j- = (1 + 7)/(l — 7) 
on the representations summed. Remark however that this equation may have no solution 
(e. g. if 7 is irrational) and should be true only in an asymptotic sense in the ultraspin 
limit where j + and j_ are both very large. 

More precisely this constraint reads 

7 >1 J± = I±i j5 n + = n_ (2.22) 

1 ±7 
7 < 1 j± = J±j = -^—^ n + = n_,. (2.23) 

where j±,j are half integers, [j 

From now on we consider only the case < 7 < 1 where the EPRL and FK models 
coincide. At 7 = 1, the EPRL/FK model reduces to a single SU(2) BF theory (see below). 



The second step replaces in each strand of (2.21) the identity lj by a projector T? 
whose definition is 



T] = d 



'u+j- fe/_/i/+=(i-7)/(i+7)] / dn\j + ,n) ® \j_,n)(j+,n\ <g> (j-,n\. (2.24) 



2 Moreover, in the case 7 > 1 the coherent states to be used below are the ones in their "anti-parallel" 
version \j,n) <g) \j,n) [9]. 

8 



Let us notice here that, in the angular momentum basis, the operator T 7 takes the form 
T ] = 5Z (U,k;j-,^j++j-,k + fy(j + +j-,m + m\j + ,m;j-,m) 

k,k,m,m 

\j+k > ®\]-~k){3+m\® < j-rh\ S m+ ^ k+h (2.25) 

where (.|.) denotes the Clebsh-Gordan coefficients. 

Grouping the four strands of a line defines a T 7 operator that acts separately and 
independently on each strand of the propagator: 

T< = ® jf ®} =1 T] f (2.26) 

so that the EPRL/FK propagator is 

C = PT 7 P; C(g, g') = f dudv f[ ]T fe_ /i/+ =(i- 7 ) /( i +7 )] a Jf (3 jf f dn f (2.27) 

J /=i if J 

Tr i /+ ®i/- (ug f (^)~ 1 t;- 1 |j /+ ,n / > ®\jf-,n f )(jf + ,n f \® < j f -,n f \) , 
where 

oij = d j+ dj_, (3j = d j++j _ (2.28) 

Lemma 2.1 The operator C is Hermitian. 
Proof We have 

Tr,-,,^, (ugjig^v- 1 ^ = fa } J dnTr jf+9jf _ (2.29) 



L 3f+$ 



ue T eg f e T e(g' f ) 1 e T ev 1 n / |j /+ j/_>0'/+J/-|nj) , 



where we have inserted the product e T e = 1 with e G SU(2) defined by (A. 8). We arrive 
at 



Tr jf+ ® jf _(ug f (g' f ) v Tj f ) = 0j f j dnTr i/+0i/ _ (we g f (g' f ) ev n f \j f+ j f _}(j f+ j f _\n}) 

= ft, J Mr^Ae'wMrfenAMt-WntoM) 
= p if j dnTr jf+9jf _ (en/ \jf+3f-){jf+Jf- \n T f e T vg' f g] l u' 1 ) 



= Tz jf+ ® jf _{vg' f g?u- l T] f ), (2.30) 

which implies the Lemma. QED 

Since the propagator is hermitian, Feynman amplitudes are again independent of the 
orientations of faces and propagators. 

Lemma 2.2 T 7 is a projector, namely (T 7 ) 2 = T 7 . 



Proof In the coherent states basis it is easier to check that (T 7 ) 3 = (T 7 ) 2 , which adding 
Hermiticity of T 7 implies the Lemma. The equation (T 7 ) 3 = (T 7 ) 2 follows from the same 
equation on each strand, since 

/3 3 / dndn'dn"\j + ,n > ®\j-,n)(j++j- } n\j++j-,n' > 



< j+ + j_ , n' \j+ + j. , n") (j + , n" | ® (j. , n" \ 

/3 2 / dndn"\ J+ , n) ® \j-,n)(j+ + j-,n\j+ + j_, n")( J+ , n"\ ® j_,n"\, (2.31) 



where we have used that 1 



3++3- 



(3,j J dn' \j+ + j- , n') (j + +j-,n'\. 



QED 



Let us also notice that the lemma is easily proven in the angular momentum basis, 



where, from (2.25) it easily follows that (T£ 



•f\2 



rr<~i 
3 ' 



Since T 7 and P do not commute, the propagator C can have non-trivial spectrum 
(with eigenvalues between and 1). Slicing the eigenvalues should allow a renormalization 
group analysis. This is why we would like to call this kind of theories dynamic GFT's. 

Remark that since T 7 is a projector, the propagator C of the EPRL/FK theory is 
bounded in norm by the propagator of the BF theory, and that Feynman amplitudes for 
the EPRL/FK theory are therefore bounded by those of the BF theory; in particular we 
expect milder ultraspin (large j) divergences in EPRL/FK. 



2.6.1 Amplitudes 

Combining the propagator and the vertex expressions, the integrations over all g, g' group 
variables can be performed explicitly, leading to the amplitude of any graph Q. This 



amplitude is given by an integral of a product over all faces of the graph as in (2.17), but 
the amplitudes for faces are different. 

To compute these face amplitudes we distinguish between closed faces (no external 
edges) and open faces (which end on external edges). 

Using the same numbering of the p edges and vertices along a closed face, its ampli- 
tude is given by 

a=1 Ua v 

where the constraint on j+,j- is implicitly understood from now on. We can perform the 
g integrals using (A.6) or (|A.7[) and we arrive at 



• / v — Z-^ a is i/+®i/-ll«=i I 



'H a ,Va"'l a+ l,V a J -3f) ' 



(2.33) 



Jf 



with hi atVa defined in (2.19) and £ p+ i = l\. Note that we use (A.6) or (A. 7) to take into 
account the fact that r)ef can change when we follow a face /. We find 






+ n ia+l,f) 



*<j/--wIC1-^L^/--w>- 



(2.34) 
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2.6.2 BF limit 



Let us see how we recover the SU(2) BF model in the limit 7 = 1. In this limit j_ = 0, 
hence j + = j + + j_ . Thus we are left with 



A 



/1 7=1 



E d T / fl *»W(J/. n iaJ \hZ:l + hl^i + \j f , n w >, (2.35) 

it a=l 



where only one SU(2) copy appears. We can use the completeness relation for coherent 



vertices and we obtain 



states (|2.14|), and cyclicity of the trace to reorder the product according to lines instead of 

(2.36) 



j J o=l 

Redefining tg a = uj + Ve a> + we finally obtain 



*4/U = E4 / ^O'nin^Cljn) = 5(II, e / 






(2.37) 



consistently with (2.8). 



2.6.3 Amplitudes with external edges 

For a face with external edges the expression is slightly modified, as there is no integration 
on the external data. 

Let us call G, G the group labels of the incoming and outgoing external strands 
respectively. We omit the edge index £ in the following. Moreover we indicate with 
M in ,fi n ,M ou t,fout the gauge transformations on the incoming, outgoing edges respectively. 
Let q be the number of internal strands. The expression of the resulting face amplitude 



is similar to (2.32) except for the fact that we don't integrate on the external labels. On 
chosing 

Vmf = Voutf = 1 (2.38) 

that is, the incoming and outgoing strand oriented according to the face, we find using 



flA.6[ ) and flA.7} 

Aext = > «,- X Tr 



E' 



u in GG l v Q ^T]u out u h T] ( Y[u ea v ea+1 TJ ) v^v-^T] 

\a=2 



(2.39) 



where, to simplify the notation, we have chosen all propagators oriented according to the 



face. It is immediately verified that it reduces to (2.33) with p = q + 2 if we glue together 
the external edges with the insertion of a delta function <5(GG _1 )5 - -/. 



11 



3 Stationary phase for BF and EPRL/FK models 

Let Q be a graph in a GFT corresponding to the BF or EPRL/FK models, made of V 
vertices, L edges and F faces, usually labelled by letters v, £ and /. In the coherent state 
basis, its amplitude can in general be written as 

i/<A J f 

where M is a normalization factor which is a rational function of the spins. As explained 
in the previous sections, the precise form of the face action and of the number of group 
variables h G SU(2) and unit vectors n G S 2 depends on the choice of the model. Note 
that the sums over the spins jf may lead to divergences, so that we introduce an ultraspin 
cut-off A that restricts the summation to spins below A. To derive the superficial power 
counting, we set jf = jkf with kf <E [0, 1] and use the stationary phase method to derive 
the large j behavior of 



[Y[dhY[dnexp{jJ2k f S f [h,n]}. (3.41) 

J f 



If the action is complex but has a negative real part, the contribution to this integral are 
quadratic fluctuations around zeroes of the real part of S which are stationary points of 
its imaginary part, otherwise the integral is exponentially suppressed as j — > oo. 

3.1 BF models in the coherent state representation 

For BF models we have one group element hi G SU(2). In dimension 4 one should work with 
SU(2) x SU(2) instead, which lead to two independent copies of the previous amplitude, 



so that we restrict ourselves to SU(2) for simplification. The amplitude is given by (2.18) 
Including the kf factor of (3.41) in the action and using 



(n,j\g\n',j) = {n\g\n') 2j , (3.42) 



with \n) a shorthand for ||n), it can be written in the form (3.40) with 



S f [h,n] = 2k f log{n\ JJ hf' f \n) (3.43) 

eedf 

Note that \n)(n\ is a projector 

\n)(n\ = l(l + (T-n), (3.44) 

so that the action reads 



S f [h,n] = fc/logTr[( J] hf f ) (l + a ■ n)] . (3.45) 



iedf 
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Since the action is the logarithm of the trace of the product of an unitary element and a 
projector, it is clear that its real part is negative (it is the logarithm of the modulus of 
the trace, obviously bounded by 1) and maximal when the unitary element is one. This is 
attained at he = 1, but other solutions may be possible. In particular, the BF amplitude 
is invariant under the gauge transformations g v at any vertex 

h e -)■ g v h t g~} (3.46) 

for any edge from the vertex v to the vertex v'. Therefore gauge transformations of the 
trivial solution he = 1 yield other equivalent solutions. More generally, there is a continuum 
of solutions connected to the trivial one which will translate into flat directions in the saddle 
point approximation. 

To perform the saddle point expansion, we expand the group element to second order 

as 

he = l-^+iA e -a + 0(A 3 i ), (3.47) 

with A e su(2) x su(2) a Lie algebra element. By the same token, we expand the unit 
vectors as 

n f =nf+^ f - fnf + 0($), with nf •£/ = (). (3.48) 

This expansion is determined by the requirement that n| = 1 up to third order terms. To 
alleviate the notation, we drop the superscript (0) in the sequel. Let us consider a face 
with edges £%,... ,£ p , then to second order 

l[h^ = l~^ + ia-A f -ia.^ f (3.49) 

leaf 

with 

A f = J2 Ve a ,fA ea and $ f = J2 Vi a jVi b j A ia A A £b . (3.50) 

l<a<p l<a<6<p 



Expanding then to second order (3.45), we get 

r Al (rifAf) 2 ■) 

S f [A £ ,Z f ] = 2k f [in f -A f --f+ y f 2 f) + i£/ • A, + in f ■ *,} (3.51) 

and we have to estimate 

/n^n^/ ex P 2 -? $>/{ in / ■ Af - ~i + {nf 2 4/) +i ^ f ■ Af +inf ■ $/ } (3 - 52) 

J £ f f 

as j — > oo. Note that we do not integrate over the vectors n/, the latter have to be chosen 
so that they are extrema of the imaginary part of S. Because all terms except the first one 
Y2f kfUf ■ Af are of second order, the imaginary part is stationary if and only if 

Y^ i k f n f -A f = ^2i Vtjkfn f -A e = MA t e M 3 , (3.53) 

/ i,f 
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which amounts to the closure condition 

Y,Vijk f n f = 0, (3.54) 

/ 

to be fulfilled for any edge £. This is the well-known requirement that, in the semi-classical 
limit, the vectors j/^/ are the sides of a triangle (resp. the area bivectors of a tetrahedron) 
that propagates along £ in dimension 3 (resp. dimension 4). The solutions of the closure 
conditions range from non-degenerate to maximally degenerate. In three dimensional (resp. 
four dimensional ) BF theory, a solution is said to be non-degenerate if all the tetrahedra 
(resp. 4-simplices) corresponding to the vertices of the graph have maximal dimension. At 
the opposite end, we say that a solution is maximally degenerate if all the vectors nj are 
proportional to a single one n , 

rif = (Tfrio with <T/ G { — 1, +1} . (3.55) 

3.1.1 Maximally degenerate case 

Let us first concentrate on the maximally degenerate solutions and show that for simply 
connected graphs (i.e. every closed loop can be shrunk to a point by deforming it through 
the faces), the quadratic saddle point approximation yields an upper bound estimate 

Ag < A 3F ~ 3r (3.56) 

with r the rank of the L x F incidence matrix rjij. This is in accordance with the general 
results for BF theory presented , see also [13] . 

To derive this result, we proceed with the following five steps. 
1. For a maximally degenerate solution, the closure constraints amount to 

£)%/*/ = (3.57) 

/ 

with Xf = kf<jf. Since the rank of the matrix r]tj is r, the Xf live in a F — r 
dimensional subspace. The signs 07 have to be adjusted so that kf > 0. We end up 



with a summation over F — r independent spins in (3.40). Let us note that since the 
incidence matrix has integer coefficients, all the spins may always be chosen to be 
half-integers, after multiplication by a suitable integer. 

2. For BF theory in the coherent state representation, we have a factor of d& per face, 
so that the normalization behaves as 

A/" = (d 3 ) 2F ~ f F , (3.58) 

where we have discarded an inessential multiplicative constant as j — > 00. 

3. The integration over £y can be performed using the Fourier representation of the 5 
function 

Jd£, exp{ij'0 • A} = j 2 5 n AA } ) (3.59) 
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with an inessential factor of (2ir) 2 absorbed in the integration measure. Note that 
the vector £f is constrained to lie in the plane orthogonal to rif, so that it enforces 
the constraint Af = only in that plane. Since Af = YleVijAi, these constraints 
are not independent. The number of independent constraints is 2r, since everything 
takes place in the plane orthogonal to a vector rif = cr/no which does not depend on 
the face. Altogether, the integration over the £/ yield a factor of j~ 2r and implement 
the constraints 

\^r]ijA£ = in the directions orthogonal to n . (3.60) 

£ 

4. Using the previous constraints, the real part of the action involving A only vanishes, 
{A f f - (n, • A f f = 0. 



5. Because the graph is simply connected, the constraints (3.60) imply the existence of 
vectors C v € R 3 attached to the vertices and orthogonal to n Q such that 

A e - (n ■ A e )n = C s{i ) - C t{e) , (3.61) 

with s(£) (resp. t(£)) the source (resp. the target) of the edge L Then, the phase 
associated to a face / reads 

n ■ $ / = ^ V£J n o ■ (C s(e) A C t(e) ). (3.62) 

ladf 

The total contribution of all the faces to the action vanishes since 

J2 k f n f -$ f = J2 Vi,fk f n f ■ (C sW A C m ) = ^ ( J2 ^f k f n f) ' { C «Q A C ^) = °» 



/■<: 



(3.63) 



using the closure condition (3.54). 



6. However, its is important to note that the components of A parallel to n are not 



constraint by (3.61) and their contribution to the action vanishes identically in the 
quadratic approximation. This is the reason why we only get an upper bound in the 
maximally degenerate case. 

Accordingly, the bound for the amplitude can be estimated as 

J2 i 2F x j~ 2r ~ A3F_3r ( 3 - 64 ) 

F — r independent spins 
of order j<A 

which is the result obtained in |12j . 

It is interesting to note that for a simply connected graph, the rank r of the incidence 



matrix can be written as r = F — (V — 1). Indeed, the system of equations (3.60), whose 
rank is 1r allows to write the 2L variables A\ in terms of 2(V — 1) differences C v — C v >, all in 
the direction orthogonal to no- Therefore, one has 2L — 2r = 2V — 2, so that r = L — V + 1 
and the amplitude of a simply connected graph scales as 

A g < A 3 ^^, (3.65) 
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with xg — F ~ L + V the Euler characteristics of the graph. This also reproduces the 
result of [13] , since xg = dimH g — dimifg + dimHg = diraHg + 1 for a simply connected 
graph. This is also in accordance with the results of [12] for graphs with planar jacket: 
The faces -F jackct of the planar jacket obey -F jacket — L + V = 2, since the associated surface 
has the topology of a sphere, and the remaining faces are in bijections with the cycles 
followed by the N = F — F jackot strands in the middle, so that the degree of divergence 
reads u g = 3(F jackot -L + ^-l) + 3iV = 3(iV + 1). 

3.1.2 Non-degenerate case 

In the non-degenerate case, the situation is slightly more complicated. The integration over 



the variables £/ yields a system of equations analogous to (3.60), but now with a vector n/ 
that varies from face to face, 

*^J rjijAt = in the direction orthogonal to Uf. (3.66) 

i 

Because for fixed £ the three (resp. four) vectors rjejnf in dimension 3 (dimension 4) span 



a space of dimension 2 (resp. 3), all the components of Ai appear in the system (3.66), 



contrary to the maximally degenerate case (3.60), which only involves the components of 



Ai orthogonal to n . In the example treated in detail below (see Sect. 3.2), (3.66) turns 
out to be equivalent to X^%/^ = 0; which has rank 3r and yields the same degree of 
divergence. In the general case, we expect non-degenerate configurations to have a less 
divergent behavior, since the degree of divergence 3F — 3r obtained in []"""] in the Abelian 
case is expected to be an overestimate in the general case and is the correct asymptotic 
behavior at least for many graphs. 

3.1.3 Two dimensional case 

To close this section, let us see how the saddle point method allows us to recover the results 
presented in [27] in the simplest case of BF theory in dimension 2. In this case GFT graphs 
are ordinary ribbon graphs with trivalent vertices representing triangles. The closure con- 
dition reads k^n^ = k/ 2 nf 2 for every edge that separates two different faces and is vacuous 
for edges that appear twice as we go along a face (if we restrict ourselves to triangulations 
of orientable surfaces). Thus, for a genus graph there is a single spin j and a single unit 
vector n. Moreover, the graph is simply connected so the phase disappears. Choosing an 
arbitrary face, all other F — 1 constraints are independent since the corresponding faces 
span a sphere with one hole that can be filled with the remaining face. Thus, r = F — 1 
and u = 3F — 3(F — 1) =3, as expected from the relation (see [2*7] ) 

Ag = J2 3 Xg ~ AXe+1 > ( 3 - 67 ) 

with xg the Euler characteristic. 

At higher genera the graph is no longer simply connected and the contribution of the 
antisymmetric part may be decisive. For instance, for the non-planar double tadpole Q\ 
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Figure 4: The "self-energy" graph Q 2 , quantum correction to the propagator. 



(torus topology) 



\-ifu \— 1> 



(3.68) 



A Gl = / dh 1 dh 2 5(h 1 h 2 (h i y 1 (h 2 ) 
an expansion to second order yields the action 

S[n,A 1 ,A 2 ]=ijn-A 1 AA 2 (3.69) 

The rank of this quadratic form is 4 and the GauBian integration over A\ and A 2 yields 

(3.70) 






f x j-" 2 ~ A, 



in accordance with (3.67). 



3.2 Divergence of the self-energy in the EPRL/FK model 

The self-energy graph Q 2 of Fig. [4] has 4 open faces. It has 6 closed faces with two edges 
each. We label the internal propagators with an index a ranging from 1 to 4 and orient 
them in the same direction. We label the 6 closed faces with pairs of indices (a, b), a < b. 
Its amplitude reads 



A 



G> 



l[dutdv a± l[A ab 



(3.71) 



a<b 



where, from (2.34) the face amplitude reads 



i 
x (j^n ah \u a -u^\j-n' ah )(j-.n' ab \v h ^v~!.\j-n ah ). (3.72) 
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\2j-| 



Using (A.17), we rewrite the amplitude above as 

3 J 

x ((n\u a -Ub±\n')(n'\v b „v-l\n)) 2j -. (3.73) 

In order to perform a stationary phase analysis we rewrite the graph amplitude as 

a ^ = E/n^n^n^n{^/)\ + ^7 ex p{^/ + +^7}} ( 3 - 74 ) 

jf a a i f 

with j r = jj^kf, kf e [0, 1] and j large. There is one coherent state per strand, which 
amounts here to label the coherent states by a couple of a face and an edge i = (/, I) such 
that r]ij y^ 0. The face action for f = ab can be written as 

Sf = 2 1 ± k f log{(n La \u^(u^- 1 \n Lb )(n Lb \v^(v^- 1 \n fta )}. (3.75) 

We employ the saddle point technique around the identity and develop the group 
elements as follows 

^ = l-^- + i^ + 0(^) 3 , ^l-^+ia-I^ + O^) 3 . (3.76) 

Morover, introducing the projector 

\ni){ui\ = -, (3.77) 

the action at the identity for the face f = ab reads 

(3.78) 

which is negative except for Uf >a = nf jb = rif. Therefore, we perform the expansion of the 
coherent state around an unit vector common to all the strands of the face 



(6 



\2 



n l = n f + Z i -^-n f + 0(£ i y, with n f • & = 0, (3.79) 

otherwise the integral is exponentially damped. To perform the expansion at second order 
of the action, it is convenient to rewrite this action as 

Sf =7 ± % log Tr{|n /i0 ) (n/, a |«a (u?) -1 \n f , b ) (n f , b \} 

+ ~f ± k f \ogTr{\n Lb )(n ffi \v b ± (v : t)- 1 \n fta )(n fta \} 
-^k f logTr{\nf, a )(n f , a \\n f , b )(n f , b \}. (3.80) 



Using the projector (3.77), the expansion to second order only involves traces of products 



of Pauli matrices and is straightforward but rather lengthy. A crucial intermediate result 
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is the expansion to second order in Ai,A 2 , £1, £2 

(£0 2 mi (^-^: 



4 



n{[l + <r(n + ft-^-n)][l- lAl ^ + i(T • (^ - A 2 + ^ A A 2 )] 
[l + cr(n + 6-^n)]} = 1 - (Al ~ ^ + in ■ (A 1 - A 2 + A, A A 2 ) 
^-^ +i(ei + e2 )^A + (ei - e2 ) ^(^-^) , (3 . 81) 



4 v^ *-/ 2 ^ ^ ' 2 

Gathering all terms together and taking the logarithm, we get 

SfiA^B^b) = k n ±{-(A±-A±f+[n r (A±-A±)} 2 (3.82) 

~(Bt ~ Bt) 2 + [n f ■ (Bt - Bt)] 2 

+in, • (A± -A} + B± - B±) +in r (A±A^ + fl? A B±) 

+i (£/,<* + e/,0 • (^ - ^ + 5± - S±) 

- faa ~^ )2 + fe,« - e w ) • [n, a k - A* - (b? - *?))] }• 

To complete the computation, one has to perform a GauBian integration with an action 
S = 2~2f(S~f + Sj)' ^ n or der to disentangle this computation, it is convenient to perform 
the following change of variables 

A± = A a ±^ Xa and B± = B a ±^ Ya . (3.83) 

The interest of this change of variables is that the terms linear in A^ 1 and B ± now only 
involve A and B, while in the quadratic terms, the pair of variables A and B on one side 
and the pair X and Y on the other side decouple. We shall return in greater detail to this 



change of variable in section 3.3 in the case of a arbitrary graph, since it allows to separate 
the action, at the level of the quadratic approximation, into a SU(2) BF action (variables 
A and B) and an ultralocal potential that only involves uncoupled variables attached to 
the vertices (variables X and Y). Turning back to the self-energy, we get 

S f {A,B,X,Y,Z] = S+lA+,B+,Z} + SjlA-,B~,$ 

= k f { - (A a - A b f + [n f ■ (A a - A b )] 2 - (B b - B a f + [n f ■ (B b - B a )] 2 
+ in r (A a -A b + B b - B a ) + in } ■ (A a A A b + B b A B a ) 
+ i fo> + £/,&) ■ {A a -A b + B b - B a ) 

~ (0 ' a ~ 2 ^ )2 + fe.a " to) " [nf A (A a -A b - (B b - B a ))] } 
+ %7 V {-{X a - X b f + [n f ■ (X a - X b )] 2 + in r (X a A X b ) 

- (Y b - Y a f + [n f ■ (Y b - Y a )} 2 + in/ ■ (Y b A y a )}. (3.84) 

Performing the GauBian integration over the two dimensional vector \f = £/,a — £/,&, one 
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has 



J d Xf expjk f { - y + X/ • [n f A (A a - A b - (B b - B a ))] } 
n } A (Aa -A b - (B b - B a )) 



2v jkf 

— exp — 

Jf 2 



(3.85) 



Discarding an inessential constant in the limit j — > oo to alleviate the notations, the graph 
amplitude can therefore be written as 

a? 2 = j2 j 18 { j n dA ° n dB - n ^ ^js BF (a, b, o 

if a a f 

]JdX a expjQ(X) x f ]JdY a exvjQ(Y) 



(3.86) 



with £/ = £/ ia + ^(,. The BF-like action is 



a<6 



71/ A (Aa -A b + B b - B a ) 

+ in ab ■ (A a -A b + B b - B a ) +in ab - (A a A A b + B b A B a ) 

+ iU-(A a -A b + B b -B a )}, (3.87) 

while the ultra local terms are 

Q[X] =- 7 + 7 - Y, Kb{ [n ab A (X a - X b )] 2 + in ab ■ (X a AX,)}. (3.88) 



a<b 



The GauBian integral over the variables A and S can be evaluated using the same tech- 



niques as in the section 34^ devoted to BF theory. Firstly, there are four closure conditions, 
one for each edge a, 



a-l 



22 k ^n a b - 22 kbaHba = ° 



(3.89) 



b=a+l 



6=1 



or explicitly, 



0. 
0. 
0. 
0. 



(3.90) 



&12™12 + k 13 n 13 + k U nu 
—kl2ni2 + ^23^23 + k 2 4n 2 4 
-^13^13 - ^23^23 + ^34^14 

^14^14 + ^24^24 + ^34^34 

Note that the first three equations are independent while the last one is their sum, so that 
the rank of the incidence matrix rjij is 3. Let us investigate the case of non-degenerate 
configurations, which means that the six vectors k ab n ab span a three dimensional space. 
Geometrically, the solution of these closure constraints can be realized by constructing a 
tetrahedron (see Fig|5]) with faces labelled 1,2,3,4 and assigning the vector k ab n ab to an 
edge between faces a and b. Consequently, we sum over 6 independent spins in (3.74). 

The GauBian integration over the variables £ ab imposes the constraints 



A a -B a 



A h -B h 



in the direction orthogonal to n, 



ab- 



(3.91) 
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Figure 5: The tetrahedron illustrating the closure condition. 



However, out of the 12 relations in (3.91), only 9 of them are independent and they are 
equivalent to 

A a -B a = C (3.92) 

with C G M. 3 . First of all, it is clear that any solution of (3.92) is a solution of (3.91). Let 



us show that the converse also holds. Let consider all the equations involving edge 1, 

A\ — B\ — A a — B a in the direction orthogonal to n la , a e {2, 3, 4} (3.93) 

Because of the closure constraint J2 a >i kia n ia = 0, the relation A\ — B\ = A a — B a holds in 
the one dimensional space orthogonal to all three vectors k\ a ri\ a . Thus the vectors A a — B a 
are all equal along this direction. We then repeat the same reasoning for the other edges 
and conclude that the vectors A a — B a are all equal along all directions using the non- 



degeneracy condition. As conclusion the rank of (3.91) is 9 since it reduces the 12 degrees 
of freedom of the 4 vectors A a — B a to a single vector and the Gaufiian integration over ^ ab 
enforcing this constraint yields a factor of j -9 . 

Using this constraint, the real part of the quadratic action obviously vanishes. The 



imaginary part can be dealt with using the techniques of section 3.1 Using first the 
constraint, we write 

= (A a + B b ) A (B a + A b ) = A a AB a + B b AA b + A a AA b + B b A B a . (3.94) 

After summation over all faces, the net contribution of the phases to the amplitude vanishes 

J2 i Khfiab ■ \K AA b + B b AB a ]=-J2 [ k abnab ■ \A a A B a + B b A A b ] , (3.95) 



a<b 



a<b 



where we have use the closure constraints. Altogether, the integration over the variables 
A, B and £ yields 



/ n dA - n dB - n ^ ^?js BF {a, b, o ~ ? 



(3.96) 
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as j — > oo. Note that this, together with a j 12 arising form the coherent state representation 
of the 5 function (j 2 per face) and a summation over 6 independent spins reproduces 

£ J 12 x r 9 ~ A 9 , (3.97) 

6 independent spins r ^ J J <CA 

which is the known result for SU(2) BF theory. Since the rank r of the incidence matrix rjij 
is 3, this reproduces with non-degenerate configurations the results of p2], with a degree 
of divergence 3F — 3r. 

Let us now consider the Gaufiian integral over the independent variables X a and Y a , 

fl[dX a expjQ(X) ~ r 1 ^ (3.98) 

J a 

which simply amounts to compute the rank of the quadratic form 

Q[X] = 7+7- Y, k ab { [n ab A (X a - X b )} 2 + in ab - (X a AX,)}. (3.99) 

a<b 

This quadratic form is associated with a symmetric bilinear form 

B[X,Z] = \(Q{X + Z,X + Z)~ Q(X -Z,X- Z)) (3.100) 

and its kernel is defined as the subspace of the variables X such that B[X, Z] = for all 
Z. First, notice that B[X,Z] is complex but the variables X and Z are real. Therefore, 
the orthogonality condition B [X, Z] = has to be fulfilled for the real and the imaginary 
part separately. Since the real part is positive (but not definite positive), X has to obey 



J2kab[n ab A(X a -X b )} 2 = 0, (3.101) 



fcafe [n a b A (A a - Ji b ) ] _ 

a<b 

or equivalently X a — X b = in the plane orthogonal n ab . Using the non- degeneracy of 
the configuration, an analysis identical to that of the constraints (3.91) leads to X a = C 
with CgI 3 that do not depend on the edge. Then, the imaginary part of the relation 
B(X,Z) = Oreads 

Y, { k ^ab ■ (C A Z b ) + k ab n ab ■ (Z a A C) } = 0, (3. 102) 

a<b 



which is identically fulfilled for any Z b because of the closure condition (3.89). Finally, the 
rank of Q is the dimension of the orthogonal of its kernel, since the latter has dimension 3 
and we have 4 vector variables X a , we obtain rank(Q) = 12 — 3 = 9, so that the Gaufiian 
integral over X yields an power of j~ 9 ' 2 . Obviously, the same holds for the integration 
over Y. 

Therefore, we obtain the power counting for the self-energy with non-degenerate 
configurations as follows 

£ j 24 xr 6 xr 9 x(r 9/2 ) 2 ~A 6 , (3.103) 

6 independent spins ~J < ^A 
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with the factor j 24 arising from a dj+dj- ~ j 2 for each of the 6 faces and a factor dj ~ j for 
each of the two strand in each face. The factor j~ 6 results form the GauBian integration 



over the 6 variables Xf — (£/> — £/,&) m (3.85) and the j from the integration over A and 



B in (3.96). This reproduces the result of [11], with non-degenerate configurations. Note 



that this is an asymptotic behavior and not a mere bound as we had before, since all the 



zero modes of the quadratic approximation correspond to gauge degrees of freedom (3.92). 

It is also of interest to notice that this result should also hold with finite non-zero 
spins on the external faces. Indeed, since the latter remain finite, the contribution of the 
external faces to the action can be neglected as j — > oo. 

Finally, let us mention that we have derived this power counting with non-degenerate 
configurations. In the next section, we shall discuss maximally degenerate configurations. 

3.3 A bound for maximally degenerate configurations 

Consider a general graph Q in the EPRL/FK model with F faces /. Since we are going 
to take the limit jf — > oo for the internal spins, the contribution of the external faces can 
be neglected, as long as their spins remain finite. Recall that the graph amplitude can be 

written as 

A g = I ]Jdh Y[A f , (3.104) 

with the face amplitude given by (2.34[). The graph amplitude may be rewritten as in 



(3.40) 



with 



A f = E { d it d ij( d i) p Jll dn U dh ^pjJ2i s i + s j}}- ( 3 - 105 ) 



Sf[n,h] = 2k n ± Yl ^{<nf*\K^* %J K^^^ lVl " 1,f M}> ( 3 - 106 ) 

1<9<P 

In the limit j — > oo, we expect that the main contribution to this integral arises from the 
neighborhood of the identity for the group elements. At the identity, the action reads 

^ ± [n,l]=2fc / 7 ± logTr[j]^(l + i^ g - ( x)]. (3.107) 



i 



This is the trace of a product of rank one projectors, its real part is always negative and 
vanishes when all the projectors are equal. Therefore, we expand the unit vectors rifj 
around an unit vector nf common to all edges of the face, 



(0, 



\2 



n f ,i q =n f + i fM - ^f^n f + 0(^ lq Y, with n f ■ £/,,, = 0, (3.108) 

together with an expansion of the group elements around the identity 

h v>l = 1 - ^- H-icr • A V>1 + 0(A Vtl ) 3 . (3.109) 
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The expansion of the action to second order follows the same steps as section 3J2 It is 
convenient to introduce 



i± 



D f.v q - e v q ,l q Vl q ,fA q ,l q + e v q ,l q+ i T Jl q +l,f A v q ,l q+1 



and 



(3.110) 
(3.111) 



®% q = Vi q jVi q+1 ,f A v q ,i q A K,i q+ i- 
After some algebra, the second order expansion of the action reads 

sf^fo] = k fl ± Yl {( n / A D i0 2 + 2[n f D t q + 2[n f ®% q 

q 
+ */* ■ (^ + ^J + */* ■ [»/ A (^ - ^J] }• ( 3 - 112 ) 



(3.113) 



In order to simplify the analysis, we perform a change of variable similar to (3.83|), 

Terms linear in A^ L are all of the form 

7 + L W) j • ^4+, + 7~L V>1 ■ A-j = L v> iA v>h with L Vjl e R 3 , (3.114) 

so that they do not involve the variables X v j Terms quadratic in A v l are all of the form 

T^^V^J+T-BK-,^], (3.115) 

where the bilinear form B[A vl , A v , v ] is either a scalar product A^-A^, v or a wedge product 
n f ' ( A vi ^ A v> /')• Substituting A^ t and A^p ;/ , it is easily seen that 

l + B[A+ l ,A+ >v ] +J-BIA-,, A- I>v ] = B[A vh A',v) +7 + fB[^,^]. (3-116) 

Then, we can express the total action Sf = Sf + Sj as a sum of a BF type action 

Sf F [A,£] = *,£ {(n, A D Ll f + 2in r D fjVq + 2in r $ fjVq 

+ i£f,i q ■ (Df^ + D f;Vq ) + £/,,, • [n, A (D f>Vq - D LVq _A] }(3.117) 



with Dr v and $/ i1; as in (3.110) and (3.111) but with A v> i instead of A vl , and an ultra 
local potential 

Q f [X\ = 7 + 7"A:/ J2 {( n f A D ^f + 2in f ■ $ /^}- ( 3 - 118 ) 
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To relate the BF face action to the more conventional one we encountered in 3.1 , let us 
perform the integration over the variables £/,„ , starting with £f tV , 

/ d£ f , Vp exp jkf{\{£f,i p ) 2 + l£f, lp ■ (0,^-x + Zf,h) + h Zfh ■ i n f A (f/4>-i - £/,0] 

£/.J, • (^/^-! + £/,«„) + £m, • [n f A (D/^ - D fiVp _ x )] } = 
I exp ^/|i(^ pi +^ /i ) + i[ nf A (e/,J,_x -^/,Ji)] 

+ i[^/^-a + £/,«, - W/ ■ (!>/,«„_! + #mJ)] + % A (D f7Vp - D f , Vp _ x )^ = 

- 2i ™/ • P/,«p_i A ^/-i) " 2(n/ A D/.^J (n f A D/.^J 

+ C/,^ • (iD f>Vp +n f A D ftVp ) + £ /jh • (!£>/,„,_! + n, A D/,^)}. (3.119) 

Note that £_f )V is orthogonal to n./ so that it couples only to the projection of Df iV t +Df jVp 
onto the subspace orthogonal to nf. Gathering all the terms in the action pertaining to 
the edges Z g _i and l±, we get 

~ 2(^,-1) _ 2~(£Wi) + 2^/.^-i ' Ovi + \ n f ' (£/A-i A €f,h) 

iO,v-i • ( D f,v P -2 + D f,v P -i + #/,<J + £/,jp_i • [n/ a p/,^ - r>/,«p_ a + £/,«„)] 

i£/,ii ■ P/,«p + D J>i + D f,v P -i) + £f,h ■ [nf A (D fm - D ftVp _ y - D fiVp j\ 
~ (D f , Vp A n f f - (D f>Vp ^ A n 7 ) 2 - 2(D f)Vp A n/) • (D fiVp _ ± A n/) 
+ 2in/ • ($„, - £>,,«,_> A D /iW1 ). (3.120) 

The integration over the variable £/,/ has a simple graphical interpretation. We have 
contracted the line l q and merged the vertex v q (associated with Djy) with the vertex v q _\ 
(associated with Df tV x ) into a new vertex (still called f g _i), associated with Df tV x +Df jVq 
and $/,«, — Df >v . A U/,^. Therefore, we may pursue this procedure till we obtain a face 



with only two edges. Then, we proceed as in section |3.2| for the self energy and integrate 
over £/i — £/,2- The remaining variable £/,i + £/,2 is a Lagrange multiplier for the constraint 
Y. q D f,v q = 0, or explicitly 

/ , e v q> i q Vi q ,f^-v q ,i q + ev q ,i q+1 Vi q +i,fAv q ,i q+1 = 0, (3.121) 



which is nothing but the constraint (3.60) written in terms of the variables e v _ lj j A v _ u i + 



e„ ^ A v t i . Then, the rest of the discussion follows the same path as in section 3.1 There 
are 2r independent constraints in the maximally degenerate case. The real part of the 
action vanishes identically as well as the imaginary part for a simply connected graph, 
once we have used these constraints and the closure constraints. Accordingly, we have 

/ 11^.' H dn f> 1 expj{ J>f [A,n]} < j<£/*/-U x r 2r , (3.122) 

J v,l f,l f 

with Lf the number of edges in the face /. 

25 



Let us finally analyze the ultra local terms given by the quadratic form Qf defined 
in (3.118[). Gathering the contributions of all faces, we get 



era = £ { I>VM"/ a h t - > r >]• + a n/ • n T a* r >}, (3.123) 

/ 

<— — > 

with Z/ ;W (resp. // )t) ) the edge entering (resp. leaving) the vertex v along the face /. First 

we notice that this is a sum over all vertices of quadratic forms defined at each vertex 

involving only variables attached to that vertex. This is the reason why we called such a 

term "ultra local". 

We then proceed as we did for the self-energy The quadratic form has a real and an 
imaginary part, but its arguments are real. Therefore, the kernel of the associated bilinear 
form is the intersection of the kernel of the real part and of the imaginary part. Because 
the real part is a sum of squares, at each vertex and for each face we have 

A «- = A -> in the direction orthogonal to n/. (3.124) 

^1 l f,v Vj ' f,v 

Since in the maximally degenerate case all the n./ are proportional to no, this simply implies 
that all vectors A v j = C v in the plane orthogonal to no, while the components collinear 



to n are left unconstrained. Then, as in section 3.2, the closure constraints imply that 
A v> i = C v also lies in the kernel of the imaginary part. If we denote by d v the valence 
of vertex v (d v may be lower than 5 since the external faces carrying spin have to be 
removed), we get a rank of 2d v — 2, (there are 3d v variables and 2 + d v solutions), so that 

jl[dXex P jQ(X) < j-^^ 1 . (3.125) 

Taking all the terms together, we get 

Ag < A 3F - 3r+F+y -£» d «. (3.126) 

The first term is the power counting of the graph Q in BF theory with group SU(2), 
while the second one results from a difference of normalization between EPRL/FK and BF 
theories. The last one is minus the half of sum of the ranks of the ultra local quadratic 
forms at each vertex. Since the latter are less important for non-degenerate configurations, 
we expect the maximally degenerate configuration to give a larger contribution, as long as 
the external spin remain finite. In particular, for the self-energy we have d v = 4 so that 
the maximally degenerate configurations are bounded by A 9 . 

When applied to the self energy with d v = 4 because we set the external spins 
to 0, we get a bound in A 9 , which suggests that degenerate configurations dominate in 
the EPRL model. However, this is only an upper bound since the zero modes of the 
degenerate configurations are not all gauge degrees of freedom, in particular the component 
of A and B along no do not contribute to the action in the quadratic approximation. 
These modes require a more thorough study involving higher order terms. Nevertheless, 
using the asymptotic behavior of 6j symbols and fusion coefficients (see [H]), we show 
in appendix [C] that degenerate configurations indeed dominate this correction to the self 
energy, but with an asymptotic behavior in A 7 instead of A 9 . This is not in contradiction 
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with the results of [TT], since the latter use the relation (6j) 2 ~ y, which implicitly assumes 
that the configuration is non-degenerate. Therefore, in the sum over spins we have to 
identify a partial sum made of spins obeying a relation such that a maximally degenerate 
configurations exist. This partial sum behaves like A 7 , while the remaining terms containing 
the non-degenerate configurations are in A 6 . 



4 Concluding remark: hint on a phase transition 

By parity, the <ft 5 4-dimensional GFT has no two point function contribution to first order 
in the coupling constant. At second order beyond the self-energy graph Q 2 , the only other 
graphs have tadpoles, hence they are absent in the colored model; in the non-colored 
model they have less faces, so we can expect the amplitude of Q2 to provide the dominant 
correction to the effective propagator of the model. 

Since that amplitude Ag 2 is positive, we can expect the whole self-energy correction 
X to be also positive. The corresponding geometric power series for the dressed or effective 
propagator 

C dressed = C + CSC + CECEC + ... = C{^^) (4.1) 

should therefore be singular when the spectrum of EC has eigenvalue 1. This should occur 
for A large enough, depending on the ultraviolet cutoff A. This is usually the signal of a 
phase transition. For instance in an ordinary 4 model a positive mass term corresponds 
to a double well potential which signals a breaking of the — > —(f) symmetry. In the vector 
<f) 4 "Ginzburg-Landau" model, it leads to the famous continuous symmetry breaking with 
appearance of an associated Goldstone boson. 

At a more speculative level, this hint of a phase transition support a scenario in 
which ordinary macroscopic smooth space-time is an emergent phenomenon. Group field 
theory, in particular its perturbative phase, might be a more fundamental description and 
space-time might result from condensation through a phase transition. This scenario is a 
version of what has been called geometrogenesis. 

In this scenario the relationship of group field theory to space-time, gravitons and 
general relativity would be somewhat similar to that between QCD and effective theories 
of nuclear forces. 
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A Harmonic analysis on SU(2) and coherent states 

We include in this appendix well known formulas for self-completeness. We start with 

dgR {j) ™(g) = 5 l0 S m0 5 k0 (A.2) 

R^™(g)RV)t (g) = £ (j,m; j, m\J,M)(J,K\j, k; ~j,k)R^{g), (A.3) 

J=h-5\ 

where R^\ (g) are unitary representations of SU(2) and (J, K\j,k;j,k) are the Clebsh- 
Gordan coefficients. We use the normalizations of [3]. 

We have 

jdgR {j) n m (g)R^%(g) = j6(jj')6?% (A.4) 

JdgR^ m n (g)R^%(g) = j-S(j,j')e^e nq (A.5) 

with i? 0) n m (#) = R^ m n (g~ l ) which imply 

J dgTr 3 Ag Tr f Bg- 1 = j6(j,f)Tr 3 AB (A.6) 

J dgTftjAgTrj.Bg = -* f 6(j,j')Tr j AeB T e T (A.7) 

with e G SU(2), 

(A.8) 



1 
-1 



We have e T e = 1 and ege T = g. 

A.l nj symbols 

We have the 3j symbols 



the 6j symbols 



(_1 )jl-h+m 3 

i mim = { ' (h,-m 3 \j 1 ,m 1 ;J2,m 2 ), (A.9) 

(Jl J2 J3 \ _ V^ „m 4 m 3 m 5 m^mx m 6 ( A 1 C\\ 

J4 J5 J6 / ^^ 

' mi..m6 

and the 9j symbols 
h h k } 

A. A A \ — V^ , m l m 2m 3 m 4 m 5 m 6 m 7 m g m 9 / A i -i \ 

J 4 Jo Jo / / » * * (, 77iim47Ti7''m2m5mg''m3m,6m9- V' rL- - L - L / 

J7 J/8 J9 J mi..m 6 

For the 15j symbols see [3]. The indices are raised and lowered with the tensor e. 
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A. 2 Coherent states 

Let us first consider the SU(2) case. We introduce the following parametrization for co- 
herent states in the spin-1/2 fundamental representation 

|I, n >=e ieA -f|-,- > (A.12) 

'2 2 2 v ' 

m = (sin0, — cos0, 0) (A. 13) 



with 



and Oi the Pauli matrices 

" 1= (! I) * 2= 'o) "=(S -l) < A " 14 ) 

This represents a rotation g n of the vector no = (0, 0, 1) of an angle 9 around the rh 
direction: 

n — > n = (sin 9 cos 0, sin 9 sin 0, cos6*) (A. 15) 

with 9 G (0, 7r),0 G (0, 27r). The coherent state \j,n > is obtained in the same way 
replacing the generators i(Ti/2 with the appropriate operators J; in the 2j + 1 dimensional 
representation. With this parametrization the scalar product of coherent states reads 

9 9 9 9 - o ■ 

(j, n\j, n) = (cos - cos - + sin - sin -e - *^ - ^) 3 (A. 16) 

which implies 

(j,n\j,n) = (\ 2 ,..| 2 

In the representation space V^' of dimension dj = 2j + 1: we have 



0>|i,fi} = (a,n|in}) 2j '. (A.17) 



li = 5^b',m)0\m|, (A. 18) 

m 

where |j, m),m G [— j,j] is the usual orthonormal basis in VK We have: 

fl^m'Gl) = O'-*"*')- ( A - 19 ) 

Hence 

C' = d, / dg R^ m {g)R {3) ^,{g). (A.20) 

JSV{2) 



B Self-energy: comparison with |11| and normaliza- 
tion conventions 

We return to the "self-energy" graph Q2 of Fig. |4| We first rewrite the propagator in a 
slightly different way, using the gauge invariance. We perform a SU(2) gauge transforma- 
tion. We multiply the u ± ,v ± variables simultaneously by SU(2) elements h, h which are 
the same for the left and right components 

u^/rV, v ± ^h~ 1 v ± (B.l) 
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and we integrate over h, h so that ( |2.27[ ) becomes: 

C(g,g') = / dhdhdudv 1 \/ j cx/Pf / dnf / 

J 3f m,rh,k,k 

y /+ t;; 1 bVm>(; + m|« + ^ + i2^ m i+ (^n / )^+^ A ((fcn/)t; 



(B-2) 



Note that we have also used (2.13). Considering the tensor product of representations 
( A. 3 ) we get 



C(g,g') 



/ dhdhdudv 1 f/^«//3/ / dnf y, (i++J-j^lJ 



+ ,m;j_,k) 



f if 



m,m,k,k 



(j+, m; j_, %+ + j_, M)RV+ + ^ M j++j _ (hn f )R^ +j -l j++j ~ ^(hn^) 



O/+Vb'+ m )0'+ m l M +#/+) ® {g'f-V-\3-k){j-k\u-gf-). 
The integration over riy through (A. 5) finally yields 



(B.3) 



C(<7,</) = J dHdndvJ[Yla Jt Tr (ug^y'v^T]^)) (B.4) 



/=i i/ 



with i? = hh>, and 

^W = &,E^W<i/+™/l«>li/-fc/> 



m. m 

kk 



< jf-k\hn f k,-M t l* f ~ kf - lSif R U ' ++Jf - )Mt M f (H) 



(B.5) 



where we have used (A. 9). Using this expression and the amplitude expression, it is checked 
below that this corresponds to the normalizations of [UJ, with k = 2. 

Note that this choice k = 2 is not the one advocated in |26j . 

To prove this statement we rewrite the graph amplitude for the "self-energy" inserting 
the new expression of the propagator (B.4). We can neglect the open faces, since the 

(B.6) 



external legs have vanishing spin. Hence, we get 

A(Q2) = / dHdudv JT Aab(u,v, H), 

J a<b 

with the face amplitude 

Aab ^2 a ja a J b d 9abdgab^j a +^j a -(u a gabgab V a lT j a ( H a)) 

L jb+ ® jb - {u b g ab ~g-^v^T] b (H b )) . 



3a,3b 

Tr, 



The amplitude for this "self-energy" graph is written in [TT] as 

a(g 2 ) = e n^*) {wt b )w; b )) 2 ([[f°) 2 > 

jab *a a<b a 



(B.7) 



(B.8) 
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where in [TT], d{j a b) — Jab i n the "ultraspin" regime. One further has 



6 J ( J 



ab) 



defined as in Eq. (A. 10), while 

j i = y' aj 12 cLj 13 cL 



Jl2 Jl3 Jl4 
J23 J24 J34 



Jl2 Jl4 Jl3 
Jl2 J 14 Jl3 
Jl2 jl4 J13 



(B.9) 



J14 



(B.10) 



and cyclically for f2, f3, fi- Both expressions, (B.6) and (B.8) are calculated for external 



j's put to zero, that is, the contributions of faces with external legs are put to 1. 



To compare our expression to (|B.8|), we perform the integration on the variables u, v 

(B.H) 



and we rewrite the integrand of face amplitudes in the form 

h 

Pab 



N ' Pab 

where we use the shorthand notation 

{g)Z = R {3)m n (9) 
We need to perform integrals of the form 



±\1ab 

mab 



%b 

m ab 



b^a 



J r, V J Vab 

b=£a 



(B.12) 



(B.13) 



(B.14) 



with a, b 



A (we have 16 integrals in total). Using (A. 3) and (A. 4) we obtain for 



dut K)» K)~ K) 



</14 



.gi29l39l4.lt 



± L,±~ l 



,±-i 



mu ± "mi2mi 3 mi4 

± _l\ fc " 



(B.15) 



dvt(vf-) (v?-) [vf-] - 4 12fcl3fcl4 ^ 2P13PM (B.16) 



_ 9 fcl2fcl3fcl4^± 
P12 ^ / P13 v 7 P14 

and similar results for a = 2,3,4. All indices are double for + and — variables, that is, 
they should carry an extra superscript (e.g. m a b — >■ m ab ). We now perform the integration 
on the variables g, g. For each face A a b they appear twice, once attached to the propagator 
containing the a variables, once to the propagator containing the b ones. This explains 
the switch in the indices below. We have 6 integrals to perform for each SU(2) copy. By 



means of (A. 4) we find 



1 



d9abd9ab{9ab9ab) kab \9ab9ab) kba = j + rf _ 

•*ab ^ab 



m ab m ba 



(B.17) 



To compare to the expression of [TT] we still have to perform the integration over H a 

(B.18) 



appearing in the T 7 . From (B.5) we have 



(n)tt(Ha) = W**M. 



Ma 



b 1a b 1a b 



(H )~ ab 

{Ua) M ab - 
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Therefore we have 4 integrals of the form 



PaY[[l P ^M r ,l„+ '"'" 



1 ab q ,q 

" ah * a 



b^a 



dH a Y\(H t 



M ab 
aJ M ah - 



b^a 



We obtain 



dH 1 (H 1 )^(H 1 )^(H 1 



)M 12 ( H )Ml,, H) M 1A = M 12 M 13 M 14 



M12M13M14 



(B.19) 



(B.20) 



and similarly for the others. We replace this result into (B.19) for each a G {1, ..,4}. We 



graph amplitude (B.6) and we obtain 



replace then all integration results (B.15), (B.16)(B.17) (B.19) into the expression for the 



Am= n^A-j ( 6 ^ 6 ^) 2 n^ 



(B.21) 



\a<b 



\a=l 



As already stated above, this reproduces (B.8) with k = 2. 



C Asymptotics of 6j and 

fusion coefficients in the degenerate case 



In this appendix, we investigate the asymptotic behavior of the self-energy correction using 
sum over spins. We first derive the asymptotics of the 6j symbols which yields the power 



counting in the SU(2) BF theory and then fusion coeffcients / appearing in (B.21) to 
obtain the power counting in the EPRL model. 



C.l Degenerate 6j and BF theory 

In the general case, the Qj symbols can be written using Racah's single sum formula (see 
for instance 1281) 



f f f) = A(l, 2, 3) A(l, 2, 4) A(l, 3, 4) A(2, 3, 4) x V ( ^ + 1)! (C.22) 
J34 Ji4 J24J ^ F(k) 



with 



A(a,b,c) 



[-jab + jbc + jacV-Uab ~ jbc + jac)Hjab + jbc ~ ja 
(jab+jbc + jac+lV- 



(C.23) 



and 



F(k) = (k-J 12 -j23-Jl3)K k -Jl3-j34-ju)K k -J23-J34-J2i)K k -Jl2-J24-juY- 

X ( J12 +323 + J34 + JU ~ k) ! (j'12 + J13 +334 +324 ~ k ) ! (j 2 3 + jl3 +324 + ju ~k)\ (C.24) 

The sums runs over all integers k such that the arguments of the factorials are non negative. 
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Consider a degenerate tetrahedron which is reduced to a single edge, whose vertices 
are labelled 1,2,3,4 with 1 and 4 on the boundary of the edge. The associated spin (lengths 
of the edges of the tetrahedron) between vertices a and b (a < b) is j a b and we have 
jac = jab + jbc if o < 6 < c. Therefore only the three spins J12, J23, J34 are independent and 
we have 

F(k) = {k- 2j 12 - 2j 23 )!(rC - 2j 12 - 2j 23 - 2j u )\{k - 2j 23 - 2j u )\(k - 2j 12 - 2j 23 - 2j u )\ 
x (2j 12 + 2.723 + 2j 34 - k)\(2j 12 + 2j 23 + 2j 3 4 - k)\(2j 12 + 4j 23 + 2j 34 - k)\ (C.25) 

The sum over k is restricted to the single term k = 2(j\ 2 + J23 + J34) and we have 

F(k) = (2j 12 )\(2j u )\(2j 23 )\. (C.26) 

There are also simplifications in the factors A, 

A (12,3)-i ^' 23 ) ! ^' 12 ) ! ^A.n.J (2j 12 )!(2j23 + 2j 3 4)! 



(2j 12 + 2j 23 + l)\ ' v ' ' ' \ (2j 12 + 2j 23 + 2j 3 4 + 1 





A(2,3,4)=l 7 ( 2 ^)'(2j 23 )! >A(1>3>4)= ( (2j34)!(2j 1 2 + 2 J 23). , ( ^ 

1 (2J34 + 2j 23 + 1)! J V ' J ^ (2j 12 + 2j 23 + 2J34 + 1)' ' 

Taking all the terms together we get 

Jjl2 J/23 jl3 1 = ,_ 1 x20' 12 +i 23 +i 3 4) 
I J14 J24J 

(2j23)!(2ji2)! (2j34)!(2j2 3 )! (2j3 4 )!(2ji 2 + 2j 23 )\ (2j 12 )\(2j 23 + 2j u ) 



X 



(2J12 + 2j 23 + 1)! (2J34 + 2j 23 + 1)! (2jl2 + 2.723 + 2j 3 4 + 1)! (2jl2 + 2j 23 + 2.734 + 1)! 



(2ji2 + 2j 23 + 2J34 + l)! (c 

(2j l2 )\(2j 23 )\(2j 34 )\ [ ■ 



which simplifies into 

ll2 323 J13 1 __ (_l)20- 12+ , 23+ ,34) 



>4 JU 324 ) V2J12 + 2J23 + V 2 j34 + 2j 23 + 1 ' 

This yields an asymptotic behavior (k a b G [0, 1] fixed) 



jh2 jk 2 3 jki3{ 1 

j^34 J&14 J^24j j^ooj 2 ' 



2 



(C.29) 



(C.30) 



Consequently, the degree of divergence of the SU{2) BF theory graph (double sunset 
without external legs) made of two vertices (each with one 6j), three edges and 6 faces 
(each with one dj = 2j + 1) is 

Regenerate = 3 + 6-2=7<d^ d ° g ° nCra t0 =9, (C.31) 

where we sum over only three spins in the maximally degenerate case. Let us note that this 
is less than the degree of divergence of non degenerate configurations, so that the latter 
are dominant in BF theory, at least for this graph. 
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C.2 Degenerate fusion coefficients and the EPRL model 

Using the notations of [UJ (Appendix B), the fusion coefficients can be expressed as a 
product of a 9j and a 3j coefficient, 



fl l (ii,j2,J3,0) = 8 i+ d + S^ tj - 8 iih y/dfr d h d j3 x 

(2j+)!(2jr)!(2j+)!(2j 2 -)^ 
(2j\+ + 2j7 + 1)!(2j 2 + + 2j 2 - + 1)!(j^^^ 

f (2j 3 + ) | (2 j 3 ~) ! Cg +j£zJL ZJL + J 2 + +J2 ) ! ^3 + +-7 3 ~ + Ji + +.7T ~ J 2 + ~ Jj ) ! \ 2 rr ooN 

^i+2i3 + +2i 3 -)!(i3 + -ii + +i2 + )Ki3 + H-ii + -i2 + )Ki3--ir+i 2 -) ' l ' ] 

with as usual j^; = —^ja- The second factor is the contribution from the 9j while the 
third one is that of the 3j. Let us notice that if 7 = 1, then j+ = j a and j~ = 0, so that 
ft % (Jii32,J3,Q) — 8 i+ ,-+5j- j-$i,j 3 - Thus, the theory reduces to a SU{2) BF theory 



.33 l ~'h 



using jB.2ip . 

The result is symmetrical in the indices 1,2,3, so let us write the degeneracy condition 
on the triangle as j\ + J2 = J3 to eliminate j 3 . After some simplifications, we get 



fi'-Uuh,h,0) = S^ }S ^ sk ^ WTWTmrTWTTy (C.33) 
Thus, the fusion coefficients scale as 

fi +i '(jki,jk 2 ,jk 3 )^ > k+ ^ O^ '* . (C.34) 

Accordingly, the power counting of the maximally degenerate configurations is (the sum- 
mation over the intertwiners i is trivial thanks to the Kronecker symbols) 

^generate = 3 + U + 3 U X (-1) + 4 X ^—\ = 7 > < p n R d L CSCnCratC = 6, (C.35) 

the first term is the contribution of the 6j (two per vertices) and the second one the 
contribution from the fusion coefficients (4 per vertices). Therefore, the degenerate con- 
figurations dominate in the EPRL model for this graph, in accordance with the quadratic 
approximation. Nevertheless, there is no reason to believe that this is a general feature of 
the model, since the quadratic approximation only yields an upper bound. 
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